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Abstract. In this work we study the phenomenon of increasing stability in 
the inverse boundary value problem for the Schrodinger equation. This prob- 
lem was previously considered by Isakov in which he discussed the phenomenon 
in different ranges of the wave number (or energy). The main contribution of 
this work is to provide a unified and easier approach to the same problem 
based on the complex geometrical optics solutions. 



1. Introduction 

Most of inverse problems are known to be severely ill-posed. This weakness 
makes it extremely difficult to design reliable reconstruction algorithms in prac- 
tice. However, in some cases, it has been observed numerically that the stability 
increases with respect to some parameter such as the wave number (or energy) (see, 
for example, [4] for the inverse obstacle scattering problem). Several rigorous jus- 
tifications of the increasing stability phenomena in different settings were obtained 
by Isakov et al [71 [H [TOl HI [2] . In particular, in [10], Isakov considered the Helmholtz 
equation with a potential 

(1.1) {A + k'^ + q{x))u{x)^OmncW 

with n > 3. He obtained stability estimates of determining q by the Dirichlet-to- 
Neumann map for different ranges of k, which demonstrate the increasing stability 
phenomena in k. The purpose of this work is to provide a more straightforward way 
to derive a similar estimate for the inverse boundary value for (|l.ip . In [10] , Isakov 
used real geometrical optics solutions for the large wave number k. In this work, 
by more careful choice of an additional large parameter and a priori constraints we 
are able to use complex geometrical optics (CGO) solutions introduced by Calderon 
[3 and Sylvester- Uhlmann fT^ for all k > 1. This will simplify the proof in |10J. 
Recently similar results were obtained by Isaev and Novikov ^ by using less explicit 
and more complicated methods of scattering theory. 
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In this work, instead of considering the Dirichlet-to-Neumann map, we define 
the boundary measurements to be the Cauchy data corresponding to ()1.1|) 

j , where u is a solution to (|1.1[) 

Hereafter, v is the unit outer normal vector of dQ. Assume that Cq-^ and Cq^ are two 
Cauchy data associated with refraction indices qi and q2, respectively. To measure 
the distance between two Cauchy data, we define 

^- .tr r \ J ■ ll(/'5) - (/,5)llffi/2eff-i/= 
dist(Cgi,CgJ = max < max mm — — , 

- (/,ff)ll_f/i/2eH-i/2 1 
max mm — — } , 

(/■9)ec,2 (/,g)ec,i ll(/,5)llffi/2e//-i/2 J 

where 

ll(/75)ll_fi-i/2e/i--i/2 = (11/11^1/2(00) + ll5llL-i/2(af2)) 
Our main theorem is stated as follows. 

Theorem 1.1. Let n > 3. Assume Cq-^ and Cq^ are Cauchy data corresponding to 
qi{x) and q2{x), respectively. Let s > n/2 and M > 0. Assume ||9;||_f/s(o) < M 
{I — 1, 2) and supp {qi — (72) C f2. Denote q the zero extension of qi — q2. Then for 
k >1 and dist(Cqj , C^j) < 1/e we have the following stability estimate: 

/ 1 \ -(2s-n) 

(1.2) Mh-hR") < Ck' dist(C,, ,Cq,) + ci^k + log ^.^^^^^^ ^ ^ 
where C > depends only on n, s, fJ, Ai^ anc? supp (gi — (72)- 



From estimate (|1.2p . it is obvious that the stability behaves more like Lipschitz 
type when k is large. We would like to point out that unlike in the acoustic case 
where the constant associated with the Lipschitz estimate grows exponentially in k 
[11] ■ the constant here grows only polynomially in k. Similarly, the corresponding 
constant obtained in [TUJ (see estimate (8) there) also grows polynomially in k. 

The paper is organized as follows. In Section 2, we will collect some known 
results about the CGO solutions and an estimate for the difference of potentials, 
which are essential tools in the proof. In Section 3, we present a detailed proof of 
Theorem 01 

2. Preliminaries 

To begin, we state the existence of CGO solutions for (|l.ip . These special so- 
lutions are first constructed by Sylvester and Uhlmann [12) . Another construction 
based on the Fourier series is given by Hahner [B] . 

Lemma 2.1. Let s > n/2. Assume that ( = rj + i^ iVy^ ^ I^") satisfies 

|77|2 = fc2 + |^|2 ?7-e = 0, 

i.e., C • C = fc^- Then there exist constants C* and C > 0, which are independent of 
k, such that if \^\ > Ci.\\q\\H^(n) then there exists a solution u to the equation (jl.ip 
of the form 

(2.1) u{x)^e'^-%l + i;{x)), 
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where ip has the estimate 



C 



Remark 2.2. Note that the correction term tJj decays in ImQ. This property is 
crucial in obtaining that the constant associated with the Lipschiz estimate grows 
only polynomially in k. 

Next inequality is an easy consequence of Alessandrini's identity. We refer to [S] 
for the proof. 

Proposition 2.3. Let ui andCq^ be solution and Cauchy data to the equation 
with q ~ qi, respectively {I — \,2). Then the following estimate holds: 



/ (92 - qi)uiU2dx 
Jn 

dui 



Ul, 



dv 



U2, 



du2 
dv 



dist(Cgi,Cgi, 



-f/i/2m_H--i/2 



3. Proof of main theorem 
To prove Theorem II. 11 we first derive two lemmas. 
Lemma 3.1. Under the assumptions in Theorem \l.l\ 



(3.1) 



\Tq{r^)\ < Ck^e'"' dist{C,„C,,) + -Mh-hr-) 



C. 



holds for k > I, r > 0, uj e M" with ^ 1 and a > C^M with P + a'^ > r^/A, 
where C > depends only on n, s, M, Q, and supp (gi — 52) and C* is the constant 
given in Lemma \2.1[ 



Proof. We will use CGO solutions (|2.ip with appropriately chosen parameter (. 
Let us denote = Vi + "i-^h I = I7 2. We can choose uj-^, u-^ G M" satisfying 

u) ■ Lu^ = UI ■ uj^ = ui^ ■ uj^ = and — \uj'^\ — 1. 



Now we set 



^1 = auj^, rji = —-UJ + \ l k'^ 
£.2 = -6 and ri2 = -ruj - rji, 

and thus 

£i-m^O: \m\' = + \£i\' 

and l^il = a > C^M > C^\\qi\\H'(n) - From Lemma [^?T1 there exist CGO solutions 

ui{x)=e'<'%l + Mx)) 
to equation p.ip with q — qi, where t/j; satisfies 

C 

Note that tpi also satisfies the estimate 

C_ 

CM CL 



(3.2) 



1^, 



i\\H^{n) 



C „ „ CM CM 
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Now, by Proposition 12.31 and using the relation —no = Ci + C21 we have that 



< 



Jn 

dist(Cqi,CqJ. 



dui 
dv 



dui 
dv 



fl'i/2e_y-i/2 



Subsequently, we obtain 



(3.3) \Tq{ru)\ 



q(a;)e-*™-^ dx 



< 



< 



g(x)e-"'^-"(l + V'i)(l+i^2)dx 

' g(a;)e-'™-^(V'i+V2 + V'iV'2)rf2: 



lii. 



dui 
dv 



Hl/20/i--l/2 



U2 



du2 
dv 



ffl/2©/i--l/2 



dist(Cg,,CgJ 



g(a:)e-™-"(^i+V2+V'iV^2)rfa; 



In view of p.3p . we want to estimate || (u/, 9it;/9j^)||^i/20^-i/2 ■ Recall that u/ 
solves (jl.ip with q ~ qi. Using assumptions ||(7z||^rs(a) < Af, and s > n/2, and 
fc > 1, we have that 

dui 



and thus 



dv 

dui 
dv 



H-i/2(an) 



<Cfc2||wi|li2(n)+C|lVui|lL2(o) 



//l/2©/i--l/2 



< Cfc2||u(||L2(S2) + C||VM;||L2(f2). 



We now choose i?o > large enough such that Vl C Bii^{Q). Then we have 
\ui{x)\ < e-«'-^(l + \ipi{x)\) < Cel«'l«« = Ce"^" 

since 

by s > n/2 and It follows that 

On the other hand, in view of \\'^'4Ji\\L^{n) < ||V''i|l-ff=(o) < C (s > n/2 > 3/2 > 1) 
and p.2p . we can estimate 

llVMnU^(o) = ||m;0 +e^^'-VV'i||i2(n) < IOIII"ilU^(o) + e'*'''^1|VV'i|U=^(n) 
< C{k + |6|)e''-"" + Cel«'l-"« = C{k + a)e''^° + Ce"^" < Cke^ 
Summing up, we obtain 



(3.4) 



ui, 



dui 
'd^ 



//l/2®ii--l/2 



<Ce\\ui\\mn)+C\\Wui\\v^-^n) 
< Ck^e^" + Cke^" < Ck 



2„Ca 



Note that here C depends on n, s, M, and the diameter of VI. 
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Let X G Co°(^) be a cut-off function satisfying x = 1 iiear supp {qi — (72), then 
we have 



g(x)e"*™'"'(V'i + V'2 + V'iV'2) dx 



(3.5) 



< / \qix)Mi'i+^2 + M2)\dx 
Jn 

< Il9llff-=(n)llx('/'i +^2 + V'l'/'2)||H=(n)■ 
Since s > n/2 and (|3.2p . we can estimate 

(3.6) ||x(V'i +i/'2 + V'iV'2)|1h=(o) 

< \\x\\H-(n){\\i'i\\H-in) + \\-ip2\\H-{n) + llV'i||ff-(n)ll'02||if=(n)) 

CM CM C CM' 
\ 1 



c 

< — . 

a 



a a a 

Finally, follows from ([331), (031), dSSJ, and ([SH). 

The following lemma is an easy corollary of Lemma 13.11 



□ 



Lemma 3.2. Suppose that the assumptions in Theorem \l.l\ hold. Let R > C^M 
with C* being the constant given in Lemma \2.1\ Then for k > \, r > Q and uj G M" 
with |w| = 1, the following estimates hold true: ifO<r<k + R then 



(3.7) 



\:Fq{rio)\ < Cfc4e^^dist(C,,,C,J + ^\\q\\H 



C. 



if r > k + R then 



(3.8) 



R' 



C. 



\Tqiruj)\ < Cfc*e^'-dist(C,,,C,J + 



Proof. It is enough to take a — R when < r < fc + i?, and take a = r when 

□ 



r > fc + i? in Lemma 13.11 

Now we prove our main theorem. 
Proof of Theorem \l.l\ Written in polar coordinates, we have that 



(3.9) 



2\-s„ri-l 



\Fq{ru)\\l + r')-'r 



diodr 



k+R 



\u\ = l 



\Tq{ruj)\\l + r^yr"--^ dujdr 



I J'g(rw)|"(l + r^)-'r"-' dojdr 



k+R J |cj|=l 



|J'q(rw)|^(l + r2)-«r"-i 



dujdr 



^■.Cih+h+h): 

where R > C^M and T > k + R are parameters which will be chosen later. 
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Our task now is to estimate each integral separately. We begin with I3. Since 
\J^q{ru!)\ < C\\qi - q2\\L^(n), qi - q2 e H^{^) and s > n/2, we get 

poo 

(3.10) ^3<Cy^ \\qi-q2\\l^n)i^ + r'r'r--Ur<CT-"^\\q,-q2\\l.^,,^ 

<CT-"(^e||Ql|2,_.(„„) + i 

for £ > 0, where m := 2s — n. 

On the other hand, by estimate p.7p . we can obtain 

(3.11) h<J^ (^CA:4e^^dist(C,,,C,J + -||g1|^-.(K„)j (H-r2)-V"-idr 



< C (^fc8e^«dist(C,,,C,j2 + _L||q^|^_^^^„^^ £ 
= C ('fc8eC^«dist(C,,,C,,)2 + l^\\q\\l_^^^ 



(l + r^)-V"-Mr 



In the same way, using estimate p.Sp . we have 

(3.12) l2<C r (ck^e^'-dist{C,„Cg,) + -\\q\\H-siR^)) (1 + r^)-^"-^ 

Jk+B. \ r J 

Jk+R 

Jk+R 



where we have used 



e'='''(l + r2)-V"-i dr < e^^ [ (1 + r^)-^"-! dr 

k+R Jk+R 
poo 

<e^* / (l + r2)-V'-idr = Ce^^, 
Jo 



T pT 



{l + r')-'r''-' dr < / r-''+''~' dr 

k+R Jk+R 



1 1 c _c 

- 2s-n + 2{k + i?)2«-"+2 - (fc + i?)2 - ;r2' 
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and s > n/2, k > I. Combining (f3J|) - (|3?T2l) gives 

(3.13) <C(/i+/2+/3) 



1 



< C (^fc»e^«dist(Q,,QJ^ + -^\\q\\H-^iM^ 

+ c(^fc«e^^dist(C,,,C,j2 + -ij||gl|^_ 

+ CfcS e'='^ dist (Cq, , J 2 + — _ — . 

To continue, we consider the following two cases: 

(i) k + R<p\og^ and (ii) k + R> p\og^, 

where R> C^M and p > are constants which will be determined later. We begin 
with the first case (i). Taking 

(3.14) i?>2VC 
and e = cT"^ (c ^ 1), we deduce that 

(3.15) Mh-^(wl'^) < Ck'^A + Ck'^e^^A + CT-^"" 

for any T > + i? by ([333]) . where A = dist{Cg,,Cg^)^. 

Now we choose T — plog(l/yl), which is greater than or equal to fc + -R by the 
condition (i) . Our current aim is to show that there exists Ci > such that 

/ -, \ -2m 

(3.16) fc^e^^A < Ci f fc + log - 



and 

(3.17) < Ci ( fc + log 



\ —27n 

A 



Substituting (|3T6| and ((3Tf| into (|3T5l) clearly implies (fr2|) . We remark that 
p.l7p is equivalent to 

(3.18) c:^/''"(^fc + logi) <plogi. 
Since we have 

fc + log < (fc + i?) + log - < (p + 1) log - 
by (i), condition p.lSp (i.e. p.l7p ) holds whenever 

(3.19) cr'/'" < 

^ ^ ^ - p+1 

On the other hand, condition p.l6p is equivalent to 

(3.20) 81ogfc+ (Cp- l)log^ + 2TOlog f fc + logi j < logCi. 
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Using (i), we can bound the left-hand side of p.20|) by 

(LHS of l[jm ) < 81ogp + 2mlog(p + 1) + {Cp - 1) log j + 2(m + 4) log log j. 
Choosing 

(3.21) p<^, 
we can see that 

(LHS of ^M)) 

< Slog ^ + 2m log + l) - ^ log ^ + 2(m + 4) log log j 

< 81og^ + 2m log ^ ^ ™>2^ (~^^ ^^"^ ''^ 4)logz 
= 8 log ^ + 2m log + 1 j + 2(m + 4) (log(4m + 16) - l) . 

Therefore, condition p.20p (i.e. p.lOp ) is satisfied provided 



(3.22) 8 log — + 2m log + 1 j + 2(m + 4) (log(4m + 16) - l) < log Ci . 

Next we consider case (ii) . We choose T — k + R and observe that the term I2 
in p. 91) does not appear in this case. Hence, instead of p.lSp . we have 



- ^ (i^ +'^"") IISIIh-mR") +^^'^''''dist(C,„C,j2 + 
Setting e = T"'/R'^ implies that 

< §IIq1I1^-=(k.) + Ck'e^^A + CR\k + Ry^^. 

Now we choose 

(3.23) R > 2VC 

and obtain that 

\\q\\jj-.^^.)<Ck^A + C{k + Rr'"^, 

which implies the desired estimate (jl.2p since from condition (ii) we have 

„ k k + R k p, 1 minjp, 1} /, , 1 

k + R> --\ > - + - bg — > ii— -!- {k + \og — 

- 2 2 -2 2^A- 2 \ ^ A 

As the last step, we choose appropriate R,p, and Ci to complete the proof. We 
first pick R > Ct:M sufficiently large satisfying p.l4p and p.23p and then choose 
p small enough satisfying p.2ip . Finally, we take Ci large enough satisfying p.l9p 
and ((3:221) . □ 
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4. Conclusion 

We think that increasing stabihty is an important feature of the inverse bound- 
ary problem for the Schrodinger potential which should lead to higher resolution 
of numerical algorithms. It is important to collect numerical evidence of this phe- 
nomenon. Our method is based on the CGO solutions constructed in jGj where the 
constants in Lemma 2.1 are explicit. So most likely one can give explicit constants 
in Theorem 1.1 at least for particular domains ft like balls. Contrary to the acoustic 
case [TT], the constants in the estimate (|1.2p depend only polynomially on k. It is 
an important and challenging question to determine whether the exponential de- 
pendence on k of the estimates in |11| is indeed generic if there are no assumptions 
on rays. 
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